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1. Introduction 

In a fully quantum-mechanical way, the spin boson model (SBM) [HO [3] is a prominent 
physical model in the research of dynamics and decoherence for numerous physical and 
chemical processes. Due to its advantage in the quantitative description of quantum 
bit (qubit), the SBM has drawn wide interest in the quantum mechanics field. In 
the last decade, many promising scalable solid-state qubit schemes have been proposed 
and realized [U El EJ [7]. Since controlling decoherence is the dominating strategy in 
solid-state qubit [H |9], qubit can be designed to be coupled to a harmonic oscillator 
(HO) or detector instead of the dissipative environment in order to minimize the 
decoherence. The HO is coupled further to the environment PEEIISEEl], which is usually 
characterized by an Ohmic spectral density Johm(^)- Such a qubit-HO-environment 
proposal can be realized as: a flux-qubit read out by a dc-SQUID [6j U\ EG] or a qubit 
placed in a leaky cavity 0, El EEE] • Since we usually only need to consider two primal 
states in the qubit and the environment is characterized by the Ohmic bath, as an 
alternative but equivalent point of view, such a qubit-HO-environment model can be 
exactly mapped to the SBM with a Lorentzian structured bath J(lj) (see ([6]) in Sec. [2]) 
[El El ESI COIl- 
Different from the Ohmic bath, the equilibrium dynamics of the SBM with such a 
structured bath or the equivalent qubit-HO-environment model is rarely studied in the 
papers, such as the studies on a driven qubit coupled to a structured environment 
by Grifoni et al. [TT1 ITS] . However, this paper is interested in the regime with 
static bias and typical studies include the quasi-adiabatic propagator path integral 
(QUAPI) [19], the Van Vleck perturbation theory together with a Born-Markov master 
equation (VVBM) [20J, the flow equation renormalization (FER) [2Tj E21 [23] , the non- 
interacting blip approximation (NIBA) [22j[2l], an approximation scheme by introducing 
a HO displacement operator [25] and the perturbation method based on unitary 
transformation [26J. It's known that QUAPI is mainly numerical and it is restricted 
to finite temperature |23j. VVBM works well with finite bias at low temperature 
and Ref. [20] is a nice and physically clear work. However, since VVBM uses the 
Van Vleck perturbation theory up to the second order in the qubit-HO coupling to 
get the eigenvalues and eigenfunctions of the non-dissipative qubit-HO system and 
solves a Born-Markov master equation for the reduced density matrix in the qubit-HO's 
eigenbasis, it requires a small qubit-HO coupling and a Born-Markov approximation 
[20J. To our knowledge, FER has not studied the non-equilibrium dynamics and it 
needs to choose extra setting parameters for its best results [27]. Until now, NIBA is 
not applicable for non-zero bias at low temperature [201 [22} [21]. Ref. [25] presents 
only results for zero bias near resonance, meanwhile, Ref. [26J with one unitary 
transformation only presents results for zero bias as well. 

In this paper, as an extension to Huang and Zheng's work J26J [28], a new 
analytical approach, beyond rotating wave approximation (RWA), based on two unitary 
transformations and the non-Markovian master equation for the density operator, is 
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applied to treat the biased spin boson model (SBM) with a Lorentzian structured bath 
for arbitrary detunings at zero temperature. One should note that the two unitary 
transformations are different from Ref. [261 128] and it makes our approach applicable 
both for non-zero bias and zero bias. Moreover, within a nontrivial Born approximation 
but without Markovian approximation, we get the analytical density operator by the 
master equation method, which can easily be extended to finite temperature comparing 
with Ref. [28J. Our approach has several advantages. First, both the localized- 
delocalized transition point oll and the coherent-incoherent transition point a c are 
studied, which have not been provided so far (except a c with zero bias by Ref. [26]). 
Second, it works well for a wide parameter range: having no direct restriction on 
the qubit-HO coupling, both for biased and unbiased, at arbitrary finite detunings 
(positive/negative detunings or on resonance) and with sufficient strong spin-bath 
coupling as long as a < a c , which is little-studied and beyond the weak coupling regime. 
Our results are checked in the exactly solvable special cases. The dynamics and the 
corresponding spectrum are compared to the literature results both for unbiased and 
biased cases. The Shiba relation and the sum rule have also been examined. 

This paper is organized as follows. Sec. [2] introduces the model and our treatment. 
Meanwhile, the ground state energy, the renormalized tunneling factor r\ and localized- 
delocalized transition point are determined. Sec. [3] presents the master equation and 
an analytical expression for the density operator without Markovian approximation. In 
Sec. HI we have calculated the non-equilibrium dynamics P(t) and the corresponding 
spectrum S(u), and presented the physical interpretation. In Sec. El it shows the 
susceptibility and the validation of the Shiba relation. The coherent-incoherent 
transition point a c is determined. 



2. Model and treatment 



In a flux qubit system, the qubit is the two macroscopically distinct quantum 
states representing clockwise and anticlockwise rotating supercurrents. The qubit is 
entangled with the detecting field, which is itself coupled with the outside noncoherent 
environment. The qubit can be characterized by Pauli matrices, the detecting 
equipment, which is actually a LC resonant circuit [61 [7] , can be described by a harmonic 
oscillator with a characteristic frequency Q and the environment can be described by a 
set of harmonic oscillators. Therefore, the qubit-HO-environment Hamiltonian can be 
written as (with Planck units h — ks — 1): 

A 



H= 



a , 



-a z + VtA^A 
2 



E 



+ (A* + A) 



go- z 



E 



(A^ + AfJ2- 



(1) 



where A is the energy difference of the qubit and e is the applied bias; A (or A*) 
and cifc (or ay are the annihilation (or creation) operators of harmonic oscillators with 
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frequencies f2 and oVs, respectively; g is the qubit-HO coupling and n k is the HO- 
environment coupling relating to the kth oscillator. The environment is described by the 
Ohmic spectral density as Johm(^) = J2k K \^{ UJ ~ UJ k) — Tcu, where Y is the dimensionless 
coupling constant to describe the Ohmic bath. 

As an alternative but equivalent point of view, such a qubit-HO-environment model 
in ([1]) can be exactly mapped to the conventional SBM with a Lorentzian structured 
bath. The Hamiltonian reads p31 M, HSl M, [26] 

H = H S + H B + H h (2) 

where the subscript 'S' denotes the spin system, the subscript 'B' denotes the boson 
environment and the subscript T denotes the interaction between the spin and the 
boson environment, with 

A e 

H B = ^u k b\b k , (4) 

k 

H i = l^Y,9k(bl + b k ), (5) 
k 

where b\ (b k ) is the creation (or annihilation) operator of the kth boson mode with 
frequency uj k \ a x and a z are Pauli matrices to describe the spin system; e is the bias, A 
is the bare tunneling and g k is the coupling between the spin and the boson environment. 
Notice that g and g k are two different quantities and there is no correlation between 
them. The boson environment of the SBM is described by the Lorentzian structured 
spectral density, and it reads [HI HH [15l [16] 

JH = ^ g k 5(u - Uh ) = {Q2 _ u2 y + {27TTujnr (6) 

in which a = lim^^o J{w)/ (2a;) = 8rcjr 2 /fi 2 is the dimensionless coupling constant. 

Our model starts with the SBM in (J2I). In order to take into account the correlation 
between the spin and bosons, a unitary transformation is applied to H to obtain H' = 
exp(S)H exp(-S), where the generator S = J2k t(W ( 2uj k))) (b[-b k )[^ k cr z + (l-£fc)cr ]. 
Here we introduce a constant cr and a fc-dependent function which will be determined 
later. We rewrite the transformed H' = H' + H[ + H' 2 as 

H' = --r]Aa x + ^ea z + ^uj k b\b k 

k 

-E^-&)-E^(i-&f, P) 

k k 

H 'i = \ E^ 1 - t>)<A + b ^ - <*) - ^^x, (8) 
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H' 2 = --Aa x (coshX - 77) - -Aia y (sinhX - rjX) 

-£|U(l-£ fe )V.-<7o), (9) 

k 

where X = J2k (gkCk/^k) (b k — bk), and 77 is the thermodynamic average of coshX, as 
r] = Z^Tr [exp(-pH) coshX] 



cxp 



^2wjf* V2T/ 



(10) 



with Z = Tr[exp(— f3H)] and T is the temperature. 

Since the spin and bosons are decoupled in H' , it is exactly solvable. By a 

unitary matrix U = (™ V \, with u = - e/W) /2, u = y/(l + e/W) /2 and 
= ^/ e 2 _|_ ^2^2^ ^he diagonalized #0 = U^H' U reads 

#0 = ~Wa z + ^u k b\b k 

k 

The eigenstate of ifo is the direct product |{rifc},±), where |±) are the eigenstates of 
o z with eigenvalues ±1 respectively, and are the eigenstates of bosons with n k 

phonons for the mode k. The ground state of H$ is |</o) — |{0fc},+) and the lowest 
exited states are ({0^}, — ), |{lfe}, +) and |{lfe}, — ). 

Similarly, we make the transformations to get H\ = WH[U and H 2 = U'H^U, 
which are treated as perturbation and they should be as small as possible. For this 
purpose, it's determined as a = —e/W and £ fc = uo k /(uj k + W). Thus 

H 1 = l(l-a z )J2Qk(bl + b k ) 

k 

+ \^2 v k [b k (o- x + ia y ) + b k (a x - ia y ) , (12) 
k 

where Q k = g k [e/(u k + W)\ and V k = g k [r}A/(u; k + W)\. In the SBM, g k < A. Q k 
and V k can be viewed as the renormalized spin-bath coupling, and they are always 
smaller than g k and even smaller for the high frequencies. Obviously, i?i|<?o) = 0. 
Under the eigenbasis of Ho, H\ has only off-diagonal terms and in the lowest states, it 
is {{%},-\H l \{l k },+) = V k , {{^ k },-\H 1 \{l k },-) = Q k and ({l fc }, -|#i|{M, +) = 
0. Meanwhile, the terms in H2 are related to the multi-boson transition and their 
contributions to the physical quantities are to the fourth order of g k (0(g%)). These are 
key points in our approach. The transformed Hamiltonian is approximated as 

H = Hq + H\ + if 2 

^H Q + H l (13) 
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in the following. Meanwhile, the previous treatment is an extension to the one proposed 
by Ref. [28J, while our generator S and the second unitary transformation are different. 
However, the fc-dependent function and decomposing the Hamiltonian into three parts 
are with the same spirit and they have been discussed detailedly in Ref. [28] . 

Thus, the ground energy of H is just the same as that of H and it is determined 

as 



The Hamiltonian H in (jSJ) can be solved exactly in two limits: one is the weak coupling 
limit with E g (a — > 0) = — \J A 2 + e 2 /2 and the other is the zero tunneling limit with 
E g (A — > 0) = — |e|/2 — Ylk9k/ (^Wk)- The ground energies in (fl4"l) are the same in both 
limits. 

Up to now, the deduction is independent of any specific spectral density and it is not 
restricted to zero temperature. In the following, the treatment is at zero temperature. 
As shown in H' Q , r\ is the renormalized tunneling factor. In the limit of zero temperature, 
it is 



The integration in (1151) can be done to the end, analytically. In the case of zero bias 
(e = 0), rj has the same expression as Huang's and similarly positive change of the 
tunneling frequency can be predicted when A ~ Q, which fails by adiabatic approach 
[26J. Generally, the renormalized tunneling factor rj is larger than 0, which means that 
there is an effective tunneling between the two states of the qubit in realistic situation. 
If the renormalized tunneling factor suddenly changes to 0, the localized-delocalized 
transition happens and the qubit will be localized in one of the two states where it is 
located before the transition. 

Fig. 1 shows numerical results of rj as a function of a. For larger T (eg. 0.15, 0.3), 
one can see that rj suddenly goes to zero at the localized-delocalized transition point 
a = (Xl, where rj = for all a > ol- While for smaller T (eg. 0.01), rj gradually goes to 
zero, and we set the cutting at the value rj = 0.0001, which is small enough. 

A phase diagram of the delocalized- localized transition point vs. bias e is 
plotted in Fig. 2 with different T (=0.01,0.02,0.05,0.1) and different detunings A/fi 
(=0.5, 1). The area of a < oil is called the "localized phase", and the area of a > oil 
the "delocalized phase". It shows that increases with increasing e and it is almost 
the same for different A's. The change of oll is remarkable for larger T or smaller e. 
Therefore, one way by applying a small bias to the qubit can be used to read out the 
already localized qubit state, since it will greatly increases and the qubit will be 
shifted from the localized state to the delocalized one. 




(14) 




(15) 
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3. Density operator and master equation 

In the Schrodinger picture, the density operator of the SBM is denoted as ps&{t) 
for the Hamiltonian H in (T5]) and the density operator for H in fflBT) is Psb(^) = 
exp(S)psB{t) exp(— S)U, where the subscript 'SB' denotes the total spin-boson 
system. In the following, it will be analyzed in the interaction picture, denoting by 
a superscript T in the operator. H is treated as the unperturbed part and Hi is really 
a good perturbed part. Moreover, in the interaction picture, it is assumed that the 
density operator for H is Pg B (t) = Ps(^)Pb, where p\{t) = Tr B Ps B (t) * s the reduced 
density operator. Within Born approximation (only keeping the second order of Hi), 
we can obtain the non-Markovian master equation for the reduced density operator 

f/s(t) = - jf T*n[ffi(t), [Hi(t'),p l s (t')p B }}dt\ (16) 

where Hi(t) is denoted as the perturbed part Hi in the interaction picture. Since the 
renormalized spin-bath coupling and in Hi are always smaller than g^, it makes 
our Born approximation nontrivial and more reasonable by comparing with the common 
used Born approximation [221 ED] , which directly does the perturbation to the second 
order of Hi in (T5]). 

The master equation in ( fl6i) . without Markovian approximation, can be done to the 
end with a Laplace transformation and an inverse-Laplace transformation. Changing 
from the interaction picture back to the Schrodinger picture, denoting the reduced 

density operator in the Schrodinger picture as ps(t) = I ^} fl ^} 2 ^\ ) for H, at zero 



P2l(t) P2 2 (t) 



temperature, we obtain 



and 



- U s _ P2i(0) f iexp{-iut)duj 

P2l[t) ~ / w -w-EM + irM' [ } 



~ (4\ — P22(0) f i exp(-ioot)duj 

P22\t) 



2tt J u- E'(lu) + ir(uj) 

— oo 

Abbreviations are used in flTjl and (JT8|) . as 



e 2 



r( w ) = 7 H + — 7 ( w -n (19) 



r] 2 A 

T] 2 A 



S(w) = R(u) + -y^R(uj - W), (20) 



and 

T'(u) = 1 (W + u) +>y(W-w), (21) 
E'(w) = R(W + u})-R(W-u), (22) 
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where R(uj) and 7(0;) are real and imaginary parts of J2k V% I { u ~ ~ u k) ( + is a 
positive infinitesimal introduced by the inverse-Laplace transformation), 

m=T ^ ] (23) 

^(uj-LUk) J (UJ - UJ )(UJ + W) 1 
k 

tf„)=,£^(„-^ = !g£W (24) 



k 



respectively. Besides, two other terms in ps{t) are p 12 (t) = [p~2i(^)]^ and pn(t) = 
1 — P~22{t)- Since the specific form of J(uj) is not involved, therefore, an analytical 
expression of the reduced density operator p~s(t) is offered and it is independent of any 
specific spectral density. 

We assume the initial density operator at t = is Psb(O) = exp(— 5*) |+) (+| pe exp(.S'). 
Thus, the corresponding initial reduced density operator for H in ffTSl is 

. l-e/W rjA/W 
,Km "I V A/W 1 + e/W I" 

4. Non-equilibrium dynamics and the physical interpretation 

For the SBM, it is common to evaluate the non-equilibrium dynamics P{t), as this is 
the quantity of interest in the experiments. P(t) is also called the population difference. 
Following the unitary transforms, it is determined as 

P(t) = Tr s (Tr B (psb(*K)) 
= Tr s (Tr B (exp (-S) Up SB (t)U^ exp (S) a z )) 

= ±[2p 22 (t) - l] + ^Re[p 21 (t)\. (26) 

Substituting (TTTl) and (TT5j) into (12"61) with the initial condition in (l2"oT) . the dynamics 
reads 

00 

2e / e \ f , cos(wt)r'(w) 

P(t) = 1 + — I / dw 9 ^ 

W ttTV V iW 7 [ u - E'(u;)] 2 + [F(a;)] 



rrW r cosM)r(^) _ ± 

^y [ w -if-sHf + [r( w )f W 







Thus, we end up at an exact analytical expression of the non-Markovian dynamics 
P(t) in (1271) . As time goes to infinity, we have the dynamics at the long time limit 
P(t -f 00) = -e/W. 
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4-1- Spectrum of the non-Markovian dynamics 

In order to get insight into the dominant frequencies of P(t), a Fourier transform is 
applied to ( 1271) according to 

POO 

S(u) = / dt cos (ut)P(t). (28) 



The spectrum S(u) is an even function and for uj > 0, it is written as 

2e / e \ r'(u) 



WJ [uj - S'H] 2 + [T'(lo)} 2 
r] 2 A 2 rV) 27re 

f ^ [ W - w - s( w )] 2 + [TM] 2 " W ^ (29) 

The frequency property of the dynamics P(t) can be analyzed by S(u) directly. 
On one hand, the first two terms in ( 1291) are Lorentzian-like functions. On the other 
hand, 7(0;) is small when g is small or u is away from Q, thus T(ui) and T'(u), which 
are functions related to j(oo), are usually small. Therefore, the dominant frequencies of 
S(u) should be the solutions u p of the equation 

u-W-E(u)=0, (30) 

and the solutions uy of the equation 

w -£'(£*;) = 0. (31) 

Since the dissipative environment generally adds shift and width to the dominant 
frequencies, we can investigate the physical nature in the limit of small HO- 
environment coupling (r — > 0). Consequently, the spectral density J{oS) in (jSJ) goes to 
{Ag 2 VL/uj) [5 (u - Q) + 5 (u + fi)] and R(u) in ([23]) goes to Ag 2 rfA 2 / [(u - + W) 2 ]. 
Therefore, according to ( 1301) and ( 1201) . the dominant frequencies u p are solutions to the 
equation 



w-W 



(n + w) 2 



n 2 A 2 e 2 

1 j (-32! 



The equation (1321) can be solved exactly. If g 2 e 2 / [A 2 (Q + W) 2 } <l,w p can be simplified 
and approximated as 



,2„2A2 



<+^_ ±V(— J (33) 

where the subscripts '1,2' relating to the sign '+, — ', respectively, and 

^V + W) + mmTw ^ m . (34) 
Similarly, according to (l3Tj) and (122]) . the dominant frequencies ay are exactly solvable, 



as 



ojp'i = (35) 
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and 

= ± f*-°f+Ww- (36) 

Since S(u) is an even function, we only need to consider the non-negative part (u > 0). 
Thus, the negative one (oy 3 ) of the solutions ui p i in fl36|) is discarded. Consequently, for 
non-zero bias, there are five dominant frequencies: u p >i = 0, u p >2, 0Jpi, uj P 2, oj p z. 
For zero bias (e = 0), we have T(u) = j(uj), = R(u) and 

r] 2 A 2 rn 

W 2 [uo-W -R(oo)} 2 + [>y(u)f 
for the non-negative part u > 0. Similarly, the dominant frequencies can be determined 
but with only two frequencies io p \ and u P 2, which are the exactly solvable solutions of 
uj-W - R(u) = 0, as 



SM = -^r ; ... , (37) 



Vt + W IfQ-WV Ag 2 7] 2 A 



,2^2 A 2 



^ = — ± i\—) + wtw (38) 

It is consistent with Huang's results (see equation (17) in Ref. [26] )• Compared to 
unbiased case, the effect of finite bias is three additional dominant frequencies: uj p >i = 0, 

UJ p ' 2, l^p3- 

Since the renormalized tunneling in the limit of small HO-environment coupling is 

2g 2 



1] = exp 



(39) 



the resonance condition is Q = Wq = A 2 + e 2 = W. For near resonance Q = W, (J33D 
can be simplified and approximated as 

TlA 

u pl , 2 = n± -L^g. (40) 

For zero bias, there are only two dominant frequencies lo p \ } 2 — according to (138|) . 
which is consistent to the result of the simple exactly solvable Jaynes-Cummings model 
[3T] . Moreover, for Q = W with finite bias, fl36|) can be simplified as 

^ 2 = ^g. (41) 

In the case of finite detunings \fl — W\ > 0, with small qubit-HO coupling 
g -C A, fl, (1331) and fl36|) can be simplified and approximated as 

4q 2 n 2 A 2 , . 

^ w+ (w-m a+ wy (42 » 

or 

n + A9 ^ 2A2 (43) 

(in (142]) and f|43|) . the larger one is uj p \ and vice versa) and 

^\w-n\ + w -iJL— . ( 44) 
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It is clear to show the physics of all these dominant frequencies: u p/1 = (13"5"|) is 
a relaxation peak, to p \ and uj p2 are related to the renormalized energy difference of the 
qubit in (jl2"j) and the energy of the HO in ( l4~3l) . and u p >2 is related to the energy difference 
of the qubit and the HO in (PHI) . Meanwhile, u p3 is related to the summation of the 
qubit energy W and the HO energy Q as shown in (1341) . Therefore, for small qubit-HO 
coupling, although u p3 and oy 2 are not exactly the summation or the difference between 
W and f2, we might still call oy 2 "beat frequency" and u p3 "sum frequency". 



Spectrum of the qubit-HO system 

Before exploring the exact spectrum S(u) corresponding to the non-Markovian P(t) in 
realistic situation, as an alternative view to (J2J), we will briefly investigate the energy 
spectrum of the equivalent qubit-HO-environment model (JT]), which is a physically 
clearer way. Since the environment generally adds shift and width to the dominant 
frequencies, to get a rough idea, the qubit-HO-environment model without HO- 
environment coupling (T = 0) is investigated here, and the qubit-HO Hamiltonian 
reads 

# q -HO = - ^o- x + e -a z + OtfA + {A ] + A)go z . (45) 

If the qubit-HO is decoupled (g = 0), (1431) is exactly solvable. By applying a unitary 
matrix to if q _Ho> h can be diagonalized as 

# q -Ho = ~<rz + VA^A. (46) 



where W = \/ e 2 + A 2 . Thus, the spectrum of the decoupled qubit-HO without 
environment is exactly shown, with eigenbasis |n, ±), where \n) denotes the eigenstates 
of HO with n (n = 0, l,---oo) phonons and |±) denotes the eigenstates of o z with 
eigenvalues ±1 respectively. 

If the qubit-HO is switched on (g ^ 0), (j45|) can be solved with exact numerical 
diagonalization, with eigenbasis denoting as \j) (j = 0, 1, • ■ ■ oo). 

To further explore the instinct of the coupled qubit-HO system, an analytical 
deduce beyond RWA is provided as follows. Since (l4"5j) has similar form to (J2]) 
when removing the summation and the multimode index k and substituting Uk — > 
Q, b k — > A\ bk — > A and gk/2 — > g, therefore, we make two similar unitary 
transformations U'' exp(S')H exp(— S')U' to H in ( 143|) with generator S' = (g/Q) (A* — 

A)[Qa z /(Q + W) - e/(Q + W')\ and U' = ( j , with v! = y/(l-e/W) /2, 

v' = a/(1 + e/W) /2 and W = \J e 2 + rj' 2 A 2 , and to the second order of the qubit-HO 
coupling g (0(g 2 )), it reaches 

ff,-Ho = - + nA ' A 
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+ ^fw' ^ A ^ ax + l(Ty) + A ^ x ~ my) ^ 

g 2 (Q + 2W) _ g 2 e 2 

(n + w) 2 n(n + w') 2, 1 ' 

where rf = exp [-2g 2 / (Q + W') 2 ] . Note that (1 - a K )\+) = 0. If the value 
\eg/ [A (Q + W')}\ < 1, then the term eg(l - a z )(A^ + A) / (Q + W) in gTD can be 
discarded. Therefore, i/q-HO is exactly solvable analytically. 

For zero bias (e = 0), Ref. [20] has used the Van Vleck perturbation up to the 
second order g and solved a Born-Markov master equation in the system's eigenbasis 
to get the dynamics P(t) = Y, n Pnn(t) + E„, ra („ >m ) Pnm(t) with the phonon number 
n, m = 0, 1, • ■ ■ oo, and it proposes selection rules for zero bias: Pnnif) vanishes 
for any n, and the non-zero p n m{t) only exists for three cases: | n even — m even | = 2, 
\n dd ~ m odd \ = 2, n even - m odd = 3 or n odd - m even = 1. When substituting e = 
into ( 1471) . the Hamiltonian ( HTt) is exactly solvable analytically. Following Ref. [20J, the 
selection rules can be deduced similarly. The selection rules show that the transition 
between the lowest energy levels |j): |0) «-> |1) and |0) «-> |2) are allowed, |1) <->• |2) and 
|0) <-> 1 3) are forbidden. This offers a second way to explain why there are only two 
dominant frequencies for zero bias. 

4-3. Results and discussion 

The Markovian approximation of P(t) is equivalent to approximate the integration in 
([IT]) and ([IB]) by the residue theorem with single pole at — 2ij and u — ij , respectively. 
It leads to 

T] 2 A 2 

P{t) = cosOot) exp (-7ot) 



47 + l)exp(-2 7o t)-l , (4* 



e 

+ W l\W 

where 7 = j(W) is the Weisskopf-Wigner approximation for the decay rate and 
ujq — W + E(W) (E(W) is the level shift). In the long time limit, the Markovian 
dynamics is the same as the non-Markovian one. 

In Fig. 3, in the case of zero bias (e = 0) with weak coupling (A = Q, g = 0A8Q, 
a = 0.004, T = 0.0154), our non-Markovian dynamics P(t) and the corresponding 
spectrum S(u) at zero temperature is compared with the ones by QUAPI [19J, by 
VVBM [20J (the numerical results) and by NIBA [21] at low temperature T = 0.1A. 
They show good agreement with both P(t) and S(u). The reasons of comparing to other 
results at the low temperature are: first, corresponding results at zero temperature by 
other methods are not found in literature; second, temperature gives a factor coth (^) 
for each frequency, and in a rough view coth (^) ~ 1 for typical frequencies (eg. to — A) 
at T = 0.1A; third, the temperature T = 0.1A is really low and the comparisons show 
that their properties are analogous. 

From Fig. 4 to Fig. 6, at finite bias (e = -0.5 A, g = 0.18 A, T = 0.0154), our non- 
Markovian dynamics P(t) and the corresponding spectrum S^u;) at zero temperature is 
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compared with the numerical results by VVBM [20] at low temperature T = 0.1 A 
for three different situations: the qubit being at positive detunings with the HO 
(ft = 1.5 A > Wo), on resonance (f2 = Wo), and negative detunings (Q = 0.9 A < Wo). 
Roughly, both the dynamics and the spectrum shows good agreement. Moreover, 
our spectrum presents five dominant frequencies: one relaxation dip (at u — 0), one 
dephasing dip (u; 2 i) and three dephasing peaks (u w , oj 2 q, oj^o), where the four damping 
oscillation frequencies are related to the energy differences of the four lowest energy levels 
of the coupled qubit-HO system as shown in the insets and they have been verified by 
the exact numerical diagonalization of -ff q _Ho (1451) . Note that the symbols lo^ denote 
the dominant frequencies of S(u) relating to the energy difference of the energy levels 
\i) and \j) of the coupled qubit-HO system. Meanwhile, the dominant frequencies can 
also be well interpreted with the relaxation dip at u> = u> p >i = in ( l35l) . the 
dephasing peaks at uio oj p2 and u 20 uj p \ in ( 1331) . uj^o uj p s in ( 1341) . and the 
dephasing dip at u 2 i 4^ u) p > 2 in (1361) . Since in Figs. 4-6 the qubit-HO coupling is small 
(g = 0.18A <C A, ft), the expressions for u p i, u) p2 and u p > 2 can be written in simpler 
approximate forms, as: for on resonance in Fig. 5, uj w <^> oj p2 and u 20 <^> oj p \ in (HO]) . 
uj 2 i ^ uj p i 2 in ( HT1) : for off-resonance in Fig. 4 and Fig. 6, uj\o & oj p2 and u 20 u p \ 
relating to the renormalized energy difference of the qubit in (1421) and relating to the 
energy of the HO in (H3|) . u 2 i <^ u p / 2 relating to the energy difference of the qubit and 
the HO in dBD. 

In Ref. [20] by VVBM, it presents four dominant frequencies: one relaxation dip 
(at uj — 0), one dephasing dip {oj 21 ) and two dephasing peaks (o^ , cu 20 ) , and similar 
result is also claimed by QUAPI in Ref. [19] (see its Fig. 7). In order to distinguish the 
dominant frequencies by different methods and/or under different conditions, analogous 
symbols uL denote the dominant frequencies of numerical results by VVBM in Ref. [20] 
are employed. As comparison, neither [19] nor [20] presents the analogous dephasing 
peak(or dip) at uj^q, and to our knowledge, it is not shown in literature. In all the three 
figures, the width and height of the dephasing peak at u 20 matches quite well with the 
one at u 20 , but our dephasing peak at uj w and dephasing dip at oj 2 \ are much higher 
and sharper, especially for on resonance in Fig. 5. Meanwhile, the dominant frequencies 
u'jj are nearly equal to the dominant frequencies Uij, but in detail uj'^ are a bit larger 
than Uij. For on resonance in Fig. 5, our dephasing dip at uj 2 \ has comparable weight 
with our dephasing peaks at c^io and u 20 , which is qualitatively different from the ones 
by VVBM. Besides, we must admit that the dephasing peaks at ^30 shown in Fig. 4 
to Fig. 6 are really small, which makes it difficult to discover. As a brief summary, 
there is a complete new dephasing peak presented in our spectrum and our dynamics 
and spectrum shows a good agreement with Ref. [20J roughly. 

In order to compare our approach with VVBM in Ref. [20], g are all rather small 
in Figs. 4 - 6, as well as the corresponding a (all a < 0.005). Nevertheless, our 
approach has no direct restriction in g and it can work with stronger a. Therefore, 
in Figs. 7 - 8, it presents with larger qubit-HO coupling (g=0.7906f2) and larger spin- 
bath coupling (a = 0.01, 0.05, 0.1) for positive detunings (Q = 2A > Wo) and near 
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resonance (Q = A ~ Wo)- The results for negative detunings are similar to the positive 
ones and the figure for negative detunings is not repeated. Likewise, the biased spectrum 
presents one relaxation peak (at uj = 0) and four dephasing peaks (c^io, ^20, ^21, ^30)- 
The dominant frequencies can be well interpreted with the energy differences of the four 
lowest energy levels of the coupled qubit-HO system as above. Similarly, they can also 
be well interpreted with the relaxation peak at uj = uyi = in fl35|) . the dephasing 
peaks at uj w <^> u p2 and uj 20 oj v \ in fl33|) . uj 30 <=> uj p3 in (1341) . and the dephasing dip 
at uj 2 i & u p '2 in ( 136|) . In contrast to the small qubit-HO coupling, the weight of the 
dephasing peaks at uj 3 o shown in Figs. 7-8 grows rather larger. For near resonance in 
Fig. 8, our dephasing peak at uj 2 i has comparable weight with our dephasing peaks at 
Uio and uj 20 and the weight of the dephasing peak at uj 30 grows much larger than that 
for off-resonance in Fig. 7. 

Meanwhile, the corresponding Markovian dynamics given by (j4"8l are presented 
in Fig. 7 to show the long time limit. Moreover, in Fig. 8 with the same A = Q, 
g = 0.7906Q and e = O.lfi, the effect with different T (= 0.002, 0.01, 0.02) is shown, 
as well as different corresponding a (= 0.01, 0.05, 0.1). The results show that the 
distributions of the dominant frequencies vary little, but with smaller r or a, the 
dephasing peaks will be higher and sharper and the dephasing will be smaller, which is 
physically reasonable. 

The non-Markovian dynamics and the spectrum for zero bias (e = 0) in Figs. 3-8 
show that the spectrum only presents two dephasing peaks (uj'{ , uj 20 ) f° r zero bias, which 
is consistent with literature results and has been interpreted in two ways as shown above 
(with (1381) or the selection rules). Similarly, u)"j are employed to denote the dominant 
frequencies for zero bias and they are related to the energy levels of the unbiased coupled 
qubit-HO system. Compared to zero bias, the effect of non-zero bias is shown in Figs. 
4-8, i.e., three more resonances in the spectrum appear: the relaxation peak at 00 — 0, 
a third dephasing peak(dip) at uj 2 \ and a fourth dephasing peak at uj 3 o- Besides, the 
dominant frequency uj" for zero bias is usually smaller than the biased one c^io and the 
dephasing peak at u'{ is usually higher. 

For non-zero bias, at uj = and LU21, it is clearly shown two dips for negative bias 
and two peaks instead for positive bias in Figs. 4-8. It can be interpreted as follows: 
uj = is mapped to oyi = and UJ21 is mapped to ccy 2 ; uyi and ujp>2 are the solutions 
ujpi of the equation (13TI) : the equation (13T|) is from the first term of S (uj) in (129]) ; the sign 
of the first term of S(uj) in ( 1291) is the same with the bias e. Similarly analysis can be 
done for the remainder three dominant frequencies. Therefore, there are always peaks 
at ujio, UJ20 and uj^o] while at uj = and u> 2 i, there are peaks for positive bias (e > 0), 
dips for negative bias (e < 0). 

As a brief summary to the biased spectrum S(uj) of the non-Markovian dynamics 
P(t): there are five resonances, i.e., the relaxation peak(dip) at u — 0, one dephasing 
peak(dip) at uj = uj 2 i, three dephasing peaks at uj = uj w , uj 2 o, ^30, which are related to 
the energy differences of the four lowest energy levels of the coupled qubit-HO system. 
For the qubit being with the HO at positive detunings (fl > Wo), the dephasing peak 
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at oj = co>io and the relaxation peak(dip) at uj = are generally dominant; for on/near 
resonance (ft = Wo), the dephasing peaks(dip) at uj = uj 10 , uj 2 o, ^21 are generally 
dominant; for negative detunings (Q < Wo), if g is small, the dephasing peak at uj = U20 
is generally dominant, otherwise, the peak at uj = uj w is dominant. A rough idea is 
that: the dominant frequency(ies) closer to the renormalized energy difference of the 
qubit (W ) usually contribute(s) more weight. 

Fig. 9 shows the dynamics with different a for non-zero bias. As usual, decay 
accompanies larger a. Our method works well for sufficient strong spin-bath coupling 
as long as a < a c (see Sec. [5]) beyond weak coupling regime. 

The sum rule of the non-Markovian dynamics is checked as shown in Table 1, and 
it is exactly satisfied with representative parameters for a < a c . 



5. Susceptibility and coherent-incoherent transition 

The susceptibility x( u ) — ~G{uj), where G{uj) (obtained in the Appendix in 
detail) is the fourier transformation of the retarded Green's function G(t) = 
— i9{t)Z~ lr Ti {exp(— (3H)[a z (t), a z }}, in which 8(t) is the unit step function. The 
imaginary part of x(^) is x"(. u ), as 

V 2 A 2 f T(lu)6(lu) 



W 2 W -E(u)] 2 + T 2 (oj) 

T(-uj)6(-u) 



(49) 



[u + W + S(-cj)] 2 + T 2 (-lu) 
and its real part x'i 00 = 0) can De obtained by Kramers- Kronig relation, as 

XV = 0) = - r ^dc^. (50) 
7T Jo w 

Our approach can be checked by the Shiba's relation [32j [33l [3H [35] 

lim^ = j[xV = 0)] 2 - (51) 

As long as a < a c , the Shiba relation is exactly satisfied as shown in Table. 1 with 
representative parameters. 

The susceptibility x"^) is the same as the second term of S{uS) in (J29"j) for oj > 
and it is an odd function of uj. Usually x"( u = 0) = 0. While increasing a to a particular 
value a c , a critical phase happens and x"{w = 0) = 00. Meanwhile, T(u) oc uj and we 
have checked that [uj — W — £(a>)] oc ^ when uj — > 0. Thus, the coherent-incoherent 
transition point [HE] cx c is defined as the solution of 

-W-E(0)=0. (52) 

In Fig. 9, non-Markovian dynamics P(t) and the susceptibility x"{ u ) are shown 
with different a (= 0.005,0.05,0.1,0.2168). Meanwhile, the coherent-incoherent 
transition point a c (= 0.21683229) is calculated by fl52|) . P(t) exhibits much abundant 
oscillation for a weak a (= 0.005), beating oscillation for a moderate a (= 0.05), a 
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badly damping oscillation for a moderately strong a (= 0.1), and nearly pure damping 
for a sufficient strong a (= 0.2168) nearly equals to a c . In the inset of Fig. 9, yf'{oS) 
is plotted against u, and the curve shows with three non-zero frequency peaks for all 
a < a c . Increasing a from weak (0.005) to strong (0.2168), the peak at the smallest 
frequency moves rapidly close to uj = and the corresponding peak grows to a great 
value, and other peaks goes close to zero. When a < a c , x"^ = 0) = 0. At a = a c , 
x"{^ = 0) = oo. Therefore, the particular value a c is the coherent-incoherent transition 
point. 

In Fig. 10, phase diagrams of the coherent-incoherent transition point a c vs. bias 
e with different T are shown. The area of a < a c is called as the "coherent phase", and 
the "incoherent phase" for a > a c . As shown in Fig. 10(a) for near- resonance A = Q, 
the changing curve of a c vs. e is an Ohmic-alike. a c gradually increases with increasing 
bias, and it remarkably increases with increasing T. As shown in Fig. 10(b) for off- 
resonance A = 0.5f2, the changing curve is nontrivial for small T = 0.075: one sharp 
peak exists around e = 0.06551) and the curve at the ends behaves Ohmic-alike. The 
sharp peak is substituted by a smooth kink for T = 0.0762. While for T = 0.08,0.1, 
the kink disappears and the whole curves behave Ohmic-alike. More results show that 
when decreasing T (eg. T < 0.075) the sharp peak grows much sharper and higher, 
and the corresponding bias of the peak becomes smaller. Under smaller T, a significant 
difference of a c between the steep area at finite bias and the platform area with zero 
bias might be utilized, eg. reading out the qubit state. 

6. Conclusions 

We have investigated the biased SBM with a Lorentzian spectral density by a new 
analytical approach at zero temperature. An equivalent description of the system is 
provided by a biased qubit coupled through a HO to an Ohmic environment. The 
starting point is the general SBM Hamiltonian (T5]) without RWA. We have applied two 
unitary transformations to the Hamiltonian and the non-Markovian master equation 
within the nontrivial Born approximation to get an expression for the density operator. 
With the density operator, we have provided analytical expressions for the non- 
Markovian dynamics P(t) and the corresponding spectrum S(uj). Meanwhile, the 
localized-delocalized transition point ai and the coherent-incoherent transition point 
a c are determined, which have not been provided so far (except a c with zero bias by 
Ref. [26]), as well as the analytical ground energy, the renormalized tunneling factor rj 
and the susceptibility x"^)- The sum rule and the Shiba relation are carefully checked 
and they are exactly satisfied as long as a < a c . 

The biased dynamics and the corresponding spectrum are key topics in this paper. 
Both for biased and unbiased, they have been compared with the results of the other 
groups and have shown good agreements. For non-zero bias, our spectrum presents five 
dominant frequencies: the relaxation peak(dip) at uj = 0, one dephasing peak(dip) at 
uj = u>2i, three dephasing peaks at uj = u>io, UJ20, UJ30, and there is a new effect: an 
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additional dephasing peak at uo = cj 30 presented in our spectrum. Our approach has no 
direct restriction on the qubit-HO coupling g. Therefore, it is a good way to investigate 
the dynamics in the little-studied strong qubit-HO coupling regime, especially in the 
static biased case, which has not been studied yet to our knowledge, as shown in 
Figs. 7-9. Moreover, the origin and the meanings of the dominant frequencies are well 
studied in two ways: providing analytical expressions for each dominant frequency in the 
limit conditions and comparing with the spectrum of the qubit-HO system. We've also 
discussed why it is sometimes peak and sometimes dip, as well as the weight distribution 
of the peaks(dips) and why there are only two dominant frequencies for unbiased. 
Meanwhile, fixing other parameters, the effect with different a and corresponding V 
is also shown. The dynamics at the long time limit is given analytically as — e/W, 
which is consistent with the Markovian dynamics. 

In summary, we have provided analytical results for interesting physical quantities 
without Markovian approximation and our approach works well at arbitrary detunings: 
on/off-resonance, with/without bias and for sufficient strong spin-bath coupling as long 
as a < a c . Admittedly, this approach is not suitable for very strong spin-bath coupling, 
e.g. a > a c . Nevertheless, the coherent regime is the most interesting one in the field 
of quantum computation and quantum information. 
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Appendix 

Following the transformation made to H to reach H, the retarded Green's function is 
G(t) = -ie{t)Z- 1 T^{exp{-(3H)[exp{iHt)a z exp{-iHt),a z }} 

= -^(t)Z- 1 Tr|exp(-^) (^[a z (t),a z ] + ^[a x (t),a x ] 

" A€ Mt),^}-^Mt),a z ]]\. (A.l) 



W 2 1 1 W 2 

The Fourier transformation of G(t) is denoted as 

e 2 2 A 2 

G H = 7J/2 a ')) + ^pT (( a x> a *)) 

-^((a z ;a x ))-^((a x ;a z )), (A.2) 

where 

((A;B)) = -ie{t)Z- x Tx {exp(-f3H)[exp(iHt)Aexp(-iHt), B] 
denotes the retarded Green's function which satisfies the following equation of motion, 
u((A;B)) = ([A,B\) + ({[A,&\;B)), 

([A,B]) = Z- 1 Tr{exp(-pH)[A,B]}. 
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Thus, we can get the following equation chain: 

u> ((cr x ; a x )) = W ((ia y ;a x )) 

k k 

u((ia y ;a x )) = 2{a z )fj Q + W ((a x ; a x )) 

+ J2® k (( a M + b k);<j x )) + Y. Vk (( a M + > ( A - 4 ) 
k k 

u {{ a A b k + h);o- x ^ = -u k ((v x (b\ - b k );a x 

+ W (iio- y {b\ + 6 fc ); a x ^ + Q k {{ia y ; a x )) , (A.5) 
u (( a x( b k ~ b k); o- x ^ = -u k ((<r x (bl + b k ); a x 

+ W (Uo- y (b\ - 6 fc ); a x )) - Q k ({a x ; a x )) , (A.6) 



+ W Ua x (b\ + 6 fc ); + Qfe (K; <r x )) , (A.7) 



w \(io- y (bl - b k ); <r x )} = -u k {(ia v (bl + b k ); a x 
+ W (U x {b{ - 6 fc ); o- x )) - Q k ((ia y ] a x )) , (A. 



where (o~ z ) ^ o = (go\o~ z \go) = 1. We have already made the cutoff approximation for the 
equation chains at the second order of g k . Besides, {{cr z ;a x )) = 0, ((c z ;a z )) = 0, and 
((o~ x ; cr z )) = 0. So the solution for G(uS) is 

^ 2 A 2 / 1 



G{u) 



W 2 \ u -w--£ k v*/(u - u k ) - Ek Qll(u -w-um) 
i 



(A.9) 
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Fig. 1 (color on the web) The renormalized tunneling factor rj vs. a. It shows 
that r) goes to zero gradually with smaller T and it goes to zero suddenly with larger T. 
A slight bias such as e = 0.0002O can be used to change the curve effectively. The inset 
is a magnification part. 

Fig. 2 (color on the web) Phase diagram of the localized-delocalized transition 
point vs. e. It shows that is almost the same with different A and the curve 
changes rapidly for larger T and/or smaller e. Parameters: near- resonance A = fi, 
r = 0.01, 0.02, 0.05, 0.1; off-resonance A = 0.5Q, T = 0.01, 0.02, 0.05, 0.1. 

Fig. 3 (color on the web) For zero bias (e = 0), non-Markovian dynamics P(t) 
and its corresponding spectrum S(uj) at zero temperature show good agreements with 
results by QUAPI [19J, VVBM [20] and NIBA [24] at low temperature T = 0.1 A. Inset 
(1) is a magnification part and inset (2) is S(uj) vs. uj. Parameters: A = Q, g = 0.180, 
a = 0.004, T = 0.0154. 

Fig. 4 (color on the web) For off-resonance Q = 1.5A with e = — 0.5A, (a) non- 
Markovian dynamics P(t) and (b) its corresponding spectrum S(u) at zero temperature 
shows good agreement with VVBM [20] at low temperature T = 0.1 A. In detail, our 
spectrum presents one relaxation dip (at lu = 0), one dephasing dip (u>2i) and three 
dephasing peaks (u w , u 2 o, co>3o)- The dominant frequency u 3Q is new. Moreover, the 
height of our peak at uiq and dip at LU21 is nearly twice the ones by VVBM. For zero 
bias, there are only two dephasing peaks. Insets: The inset in (a) and insets (1)~(3) in 
(b) present magnification parts. Inset (4) in (b) presents the schematic energy levels of 
the qubit-HO system with finite bias, where the dashed lines on the left part for the 
uncoupled (g = 0) and the solid lines on the right part are numerical exact results for the 
coupled. Our dominant frequencies are related to the energy differences of the coupled 
qubit-HO system with numerical calculation. Parameters: g = 0.18A, T = 0.0154, 
a=0.00177. 

Fig. 5 (color on the web) For on-resonance f2 = Wq = VA 2 + e 2 with e = — 0.5A, 
(a) non-Markovian dynamics P(t) and (b) its corresponding spectrum S(u) at zero 
temperature shows good agreement with VVBM [2U] at low temperature T = 0.1 A. 
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Similarly to Fig. 4, our spectrum presents five dominant frequencies (oj = 0, u 10 , u 20 , 
^21, ^30) an d the dominant frequency u 30 is new. Our dephasing dip at u 2 i is much 
deeper and sharper than the ones by VVBM. Our dip at uj 2 \ has comparable weight 
with our dephasing peaks at Ui and u 20 , which is qualitatively different from the ones 
by VVBM. For zero bias, there are only two dephasing peaks. Insets: The inset in (a) 
and inset (1) in (b) presents magnification parts. Inset (2) in (b) presents the schematic 
energy levels of the qubit-HO system with finite bias, where the left dashed lines are for 
the uncoupled and the right solid lines are for the coupled. Our dominant frequencies 
are related to the energy differences of the coupled qubit-HO system with numerical 
calculation. Parameters: g = 0.18A, T = 0.0154, a=0.00319. 

Fig. 6 (color on the web) For off-resonance Q = 0.9A with e = — 0.5A, (a) non- 
Markovian dynamics P(t) and (b) its corresponding spectrum S(u) at zero temperature 
shows good agreement with VVBM [20] at low temperature T = 0.1A. Similarly to 
Fig. 4, our spectrum presents five dominant frequencies (a; = 0, uj W) u 20 , uj 2 %, w 30 ) and 
the dominant frequency uj^q is new. Moreover, our peaks (or dips) are a bit sharper 
and higher than the ones by VVBM. For zero bias, there are only two dephasing peaks. 
Insets: The inset in (a) and inset (1) in (b) presents magnification parts. Inset (2) in 
(b) presents the schematic energy levels of the qubit-HO system with finite bias, where 
the left dashed lines are for the uncoupled and the right solid lines are for the coupled. 
Our dominant frequencies are related to the energy differences of the coupled qubit-HO 
system with numerical calculation. Parameters: g = 0.18A, T = 0.0154, a=0. 00493. 

Fig. 7 (color on the web) For off- resonance A = 0.5Q with larger qubit- 
HO coupling ^=0.790647, (a) non-Markovian dynamics P(t) and (b) its corresponding 
spectrum S(u) at zero temperature is shown with bias e = 0.1Q and with zero bias, 
while the Markovian ones in (a) has shown the long time limit. Similarly, the biased 
spectrum presents four dephasing peaks (uj w , u 20 , u> 2 i, ^30) an d one relaxation peak 
(00 = 0), while the unbiased one only presents two dephasing peaks. Insets: Inset (1) in 
(b) presents a magnification part. Inset (2) in (b) presents the schematic energy levels of 
the qubit-HO system with finite bias, where the left dashed lines are for the uncoupled 
and the right solid lines are for the coupled. Our dominant frequencies are related 
to the energy differences of the coupled qubit-HO system with numerical calculation. 
Parameters: a = 0.01, T = 0.002. 
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Fig. 8 (color on the web) For near-resonance A = VL with larger qubit-HO coupling 
c/=0.7906Q, (a) non-Markovian dynamics P(t) and (b) its corresponding spectrum S(u) 
at zero temperature is compared between the biased e = 0.1Q and the unbiased at 
a = 0.05, T = 0.01. Similarly, the biased spectrum presents four dephasing peaks (a>io, 
^20) ^21, ^30) an d one relaxation peak {uj — 0), while the unbiased one only presents two 
dephasing peaks. Besides, with the same A = Q, g = 0.7906f2 and e = 0.1Q, the effect 
with different r or a (red dot line: a = 0.01, T = 0.002 and violet dash line: a — 0.1, 
T = 0.02) is shown: the distributions of the dominant frequencies vary little, but with 
smaller T or a, the dephasing peaks will be higher and sharper and the dephasing will be 
smaller. Insets: Inset in (b) presents the schematic energy levels of the qubit-HO system 
with finite bias where the left dashed lines are for the uncoupled and the right solid lines 
are for the coupled. Our dominant frequencies are related to the energy differences of 
the coupled qubit-HO system with numerical calculation. 

Fig. 9 (color on the web) Non-Markovian dynamics P(t) for the near resonance 
case A = f2 with e = O.lfi and T = 0.01 for different spin-bath coupling a = 0.005, 
0.05, 0.1, 0.2168 (corresponding qubit-HO coupling g/Q = 0.25, 0.7906, 1.1180, 1.6462) 
with all a < a c = 0.21683229. It shows that when increasing a from weak to strong, 
the dynamics goes from abundant oscillation to nearly pure damping. Insets: The inset 
is the susceptibility x"^) vs - 00 an d it shows that when increasing a close to a c the 
highest peak goes close to infinity near uj = 0. 

Fig. 10 Phase diagram of the coherent-incoherent transition point a c vs. bias 
e, for (a) near-resonance A = with different T = 0.05, 0.08, 0.1; (b) off-resonance 
A = 0.5fi with different T = 0.075, 0.0762, 0.08, 0.1. A kink appears when T is small 
for off-resonance, otherwise a c gradually increases with e. 
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Table 1 The sum rule and the Shiba relation are checked with representative 
parameters. Here R = lim^o F(uj)/j [x'(0)] 2 , where F{uj) = x"(u;)/J(u;). 
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